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'. Abstract 

& | We consider 5-dimensional cosmological solutions of a single brane. The correct cos- 

mology on the brane, i.e., governed by the standard 4-dimensional Friedmann equation, 
and stable compactification of the extra dimension is guaranteed by the existence of a non- 
vanishing T| which is proportional to the 4-dimensional trace of the energy-momentum 
tensor. We show that this component of the energy-momentum tensor arises from the 
backreaction of the dilaton coupling to the brane. The same positive features are exhib- 
ited in solutions found in the presence of non-vanishing cosmological constants both on the 
brane (A^.) and in the bulk (A#). Moreover, the restoration of the Friedmann equation, 
with the correct sign, takes place for both signs of A# so long as the sign of Ab r is opposite 
A^ in order to cancel the energy densities of the two cosmological constants. We further 
extend our single-brane thin-wall solution to allow a brane with finite thickness. 



1 Introduction 



It goes without saying that there has been an exceptional amount of interest in theories 
with extra dimensions, and in particular the resulting cosmologies that arise in those the- 
ories. In the past, in theories with extra dimensions, one often simply assumed a small 
compact extra dimension, of order the Planck size. In this case, it is usually sufficient to 
average over the extra dimension(s) and once fixed, cosmology on the large 4D space-time 
is conventional, other than the appearance of additional massive degrees of freedom. In 
10D string theory, while the structure of the compact 6D space has a profound effect on 
the resulting interactions of matter in the effective low-energy theory, cosmology is only 
affected by terms suppressed by inverse powers of the Planck scale. Below the Planck scale, 
these are generally negligible. When string theory is promoted to M-theory, the size of 
the extra 11th dimension, is generally assumed to be somewhat larger than Mp 1 , and thus 
allows for a separation of the string (unification) scale and the Planck scale |jj . 

By allowing the extra dimension to be significantly larger than the Planck scale, one 
can try to relate the electroweak scale to the fundamental higher dimensional Planck scale 
P, |3|, H|]. Alternatively, it has been suggested that an exponential scaling of the metric 
along the extra dimension (or "warp" factor) can also lead to a resolution of the hierarchy 
problem ||. In both of these approaches the extra dimension may be as large as a mm 
scale (in the former) and even infinite (in the latter, see also ||). In these theories, it is 
necessary to assume that standard model particles are confined to a 3+1-dimensional brane 
of an n+3+l-dimensional space-time. On the other hand, gravity propagates in the bulk. 

Many theories have been constructed using the notion that there are in fact two branes 
which represent the boundaries of a higher dimensional space-time. Such was the motivation 
from the Horava- Witten model in M-theory |1| . Indeed, there has been considerable effort 
to produce cosmological models on such boundaries |7], ||, [|. In a different set of solutions 
0, the brane is a slice of AdS and the negative cosmological constant of the bulk is used to 
cancel the cosmological constant or tension on the brane. In the single brane scenario |J, 
gravity is effectively confined to the brane by the steep warp factor generated by the tension 
dominating the brane. In most other constructions, two branes are necessary to achieve 
compactification. Of course in realistic cosmologies, the energy density of the Universe 
must be dominated by matter. 

An important observation about the cosmology of the "brane world" was made in Refs. 
H, [T(|. These papers showed that the Hubble parameter H governing the expansion 
of the scale factor on the brane has a different behavior than derived from the usual 4-d 
Friedmann equations. In particular, the Hubble parameter is proportional to the energy 
density on the brane instead of the familiar dependence H ~ ^fp. While mathematically 
correct as a solution to Einstein's equations, this behavior, if nothing else, indicates that 
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some key ingredient is missing if this theory is to be capable of describing our Universe. It 
is worth noting that this unusual behavior of H does not depend on the size of the extra 
dimension, but holds even for r ~ M~ l ~ Mp^. While this abnormal behavior of H may 
be possible at very early cosmological epochs which cannot be immediately confronted with 
observational data, the H ~ p -type of behavior is totally unacceptable at later times, 
particularly at big bang nucleosynthesis. 

The origin of the aberrant expansion law can be traced to a condition, imposed on the 
solution of the 5D Einstein's equations, which relates the energy density on the two branes. 
Roughly speaking, the presence of matter on the 3-brane, produces a y-dependence (y is 
the spatial coordinate of the extra dimension) of the 3-space scale factor a such that a(y) 
decreases as one moves away from the brane. Indeed for an infinitely thin brane, there is 
a discontinuity in the derivative of a with respect to y. (This discontinuity is removed by 
considering branes of finite thickness, but the qualitative behavior of the solutions are un- 
changed as we will show below.) For a brane dominated by a positive tension, a(y) decreases 
exponentially away from the brane, and in effect compactification is unnecessary ||. If the 
brane is not dominated by the tension, the decrease in a(y) is softer and compactification 
may be necessary. To compactify, a second brane must be placed at some fixed distance 
from the initial brane. The scale factor will also experience a discontinuity (in a'(y)) on the 
second brane, but of opposite sign. This distinction imposes very severe constraints on the 
matter and tension of the second brane. In fact, if the energy density on the original brane 
is pi, then on the second brane it must be p 2 = —p\ + 0(p 2 ). At cosmological distance 
scales (much larger than the separation of the branes), the Hubble parameter will be given 



by H Mp — pi + p 2 ~ p . This point was recently emphasized in [] I 



A number of "remedies" to this problem have been suggested [12], [L3|, [14], [L5| . Most of 
these works exploit the Randall- Sundrum construction with two branes of opposite tensions 
±A& r and obtain the correct cosmological Friedmann equation by cancelling the leading A^ r - 
proportional term with a negative bulk cosmological constant A#. The remaining cross-term 
leads to the normal expansion law and the 4D Plank scale is related to the tension on the 
brane. For an alternative solution to this problem see [HJ. Cosmology of the single-brane 
scenario of ref . has been considered in JT3], [H| . For additional papers of interest related 
to the Randall-Sundrum scenario see [JTSj . 

In a different approach to the problem |H| we derived sufficient conditions which ensure 
a smooth transition to conventional cosmology and Newton's law on the brane. We found 
that if Tj^ lk ^ 0, and in particular, T 55 ^ on the brane and in the bulk, the transition to 
the conventional cosmology can be obtained (for both thin and thick branes). Furthermore, 
the same class of solutions [[HJ allowed for the possibility of compactification with single 
brane, and obviates the need for an exotic object such as a brane with negative matter 
density. 
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The main result achieved in |19| is that the desired resulting expansion law (H 2 oc p), 
together with a fixed size (L = const) for the extra dimension can be achieved when T55 7^ 0. 
In fact, the specific value of T| must be chosen to be proportional to the trace of the 4D 
energy-momentum tensor on the brane and inversely proportional to the size of an extra 
dimension, T| oc (— p + 3p)/2L + 0(p 2 ). It was hypothesized in |L9j that this value is in 
fact the consequence of the physics responsible for the stabilization of the dilaton. 

In this paper we take this idea one step further and study the relation between T| 
and the dilaton fixing in more detail. In the absence of the dilaton stabilization mecha- 
nism certain conditions should be imposed on the 5D energy-momentum tensor to preserve 
dilaton flatness and thus the stability of the radius of an extra dimension. Working to 
linear order in energy density (assumed small), we find a generic integral condition on the 
energy-momentum tensor of 5D matter. When the explicit dilaton potential is introduced, 
we show that T| oc (— p + 3p)/2L + 0(p 2 ) arises as a back-reaction of the dilaton potential 
in the presence of the brane and, to this order, does not depend on the details of the dilaton 
fixing potential. That is, the solution found in [19] is completely generic. 

Another important goal that we pursue in this paper is the development of our initial 
solution [19] in several new directions. In particular, we study whether the thin wall solution 
with a compact extra dimension allows for the generalization on the case of small but finite 
thickness. We also extend our solution to include cosmological constants on the brane and 
in the bulk, without restricting the sign of the bulk cosmological constant. In all of these 
cases we recover normal Friedmann-type solutions on the brane by requiring the stability 
of the transverse dimension. 

We organize this paper as follows. The main ingredients of the model, Einstein equations 
specialized for particular ansatzes for the metric and energy-momentum tensor, are intro- 
duced in the next section. Section 3 deals with the interpretation of the (55)-component 
of the energy- momentum tensor in the bulk. First, we consider the condition for dilaton 
flatness in the presence of 5D matter and show the desirable form of T| which satisfies this 
condition. Then we turn to the most important case when the dilaton receives an explicit 
potential and demonstrate how T| is generated as the response of this potential to the pres- 
ence of the brane. In section 4, the thin wall solution is generalized to take into account 
cosmological constants on the brane and in the bulk. We address the cases of negative, 
positive and zero cosmological constant in the bulk. In section 5 the single brane solution 
with compactification is analyzed when the brane is given a certain thickness A C L. Our 
conclusions are presented in section 6. 
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2 The Theoretical Framework 



As an example of a higher-dimensional theory describing the coupling of the matter content 
of the universe with gravity, we consider the following 5-dimensional theory 



S = Jd 5 x^G^{^R + C } 



(2.1) 



where C a represents all possible contributions to the action which are not strictly gravi- 
tational. In the above, M5 is the fundamental 5-dimensional Planck mass, and the hat will 
denote 5-dimensional quantities. The line-element of the 5-dimensional manifold is given 
by the following ansatz 



ds 2 



-n 2 (t, y)dt 2 + a 2 (if:, y)8 i jdx l dx^ + b 2 (t, y)dy 2 



(2.2) 



where {t,x 1 } and y denote the usual, 4-dimensional spacetime and the extra dimension, 
respectively. 

When the action functional ( |2.1| ) is varied with respect to the 5-dimensional metric 

(5) 

tensor G MN , Einstein's equations are derived, which for the above spacetime background 
take the form (see e.g. 
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where k 2 = 8nG = 87r/M| and the dots and primes denote differentiation with respect to 
t and y, respectively. 

The 5-dimensional energy-momentum tensor of the theory is obtained by the variation 
of y—G^> L with respect to the metric and may be written in the form 



r M ^ = diag(-p,p,p,i),T 5 & ) 



(2.7) 
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and, when combined with the equation for the conservation of the energy-momentum tensor, 
DmT m n = 0, leads to the following relations 




(2.9) 



(2.8) 



In what follows, we adopt the brane-world scenario and, thus, we assume that the usual 
matter content of the universe is confined on a brane which may have zero or non-zero 
thickness. However, we make no such assumption regarding the fifth component T| which 
can be smoothly distributed along the extra dimension while non-vanishing p and p may 
exist in the bulk due to a cosmo logical constant. 

In the framework of the above theory, we are looking for cosmological solutions which 
admit compactification of the extra dimension and, at the same time, respect the usual 
form of the 4-dimensional Friedmann equation that arises from Einstein's equations in 
the limit y — > 0. Since we are interested in a compact extra dimension, we impose the 
condition that the scale factor of the fifth dimension depends neither on space nor time, 
i.e. b = bo = const. In that case, our task is greatly simplified since the (05)-component 
(|2.5| ) of Einstein's equations can be easily integrated to give the result 



Moreover, the (OO)-component ( |2.3| ), now, reduces to a second-order differential equation 
for a(t,y) with respect to y with the general solution depending on the form of the energy 
density p of the universe. The function X(t) is of equal importance as it leads, not only to 
the determination of the lapse function n(t,y), but also to the 4-dimensional Friedmann 
equation at the origin, that is, on the brane. Using the normalization n(t,y = 0) = 1, the 
Hubble parameter can be expressed in terms of X(t) in the following way 



where the subscript denotes quantities evaluated at y = 0. As it will be demonstrated 
shortly, the form of the function A(t), and thus the restoration of the Friedmann equation 
at the origin, will strongly depend on the expression of the fifth component of the energy- 



It is important to remark at this point, that so far, we have considered the constancy 
of the scale factor b as an external condition. That is, it is not one which is obtained 
automatically from the solution of Einstein's equations. Fixing b, is equivalent to fixing 



n(t,y) = X(t) a(t,y) . 



(2.10) 





momentum tensor T|. 
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the vev of the dilaton in the dimensionally reduced theory. This is clearly a physically 
reasonable assumption since a rolling dilaton implies the variation of gauge couplings and 
particle masses and there are very strong constraints against this p0[ . Indeed one would 
expect that a dilaton potential is generated elsewhere in the theory (buried in C Q ) and 
provides for the necessary additional condition which fixes b. As we shall see in the next 
section, such a potential automatically gives rise to the T 55 component advocated in fl9fl. 



3 Stabilization of the extra dimension 

Suppose for a moment that we consider the simple example of a 5D Kaluza-Klein theory 
without matter. As such, we will have a solution for any size of the compactified dimension, 
i.e. the dilaton potential is flat. Let us now suppose that 5D matter with a given energy- 
momentum tensor is distributed in this space. Clearly, there is no guarantee that a constant 
radius will remain an integral of motion unless the energy-momentum of the matter satisfies 
certain constraints. We next consider the conditions which must be imposed on T M N in order 
to preserve the dilaton flatness. 

In some sense the answer to this question is similar to the following toy example. Imagine 
a massless scalar field 4>(x) in d dimensions which interacts with an external current J(x). 
The condition 

' J(x)d d x = 0. (3.1) 



ensures the stability of any vacuum configuration (<fi(x)) — v at the classical level. An 
analogous integral condition exists in the case of the Kaluza-Klein theory. The dilaton 
remains in equilibrium and the size of an extra dimension stays fixed if the stress energy of 
the 5D matter satisfies the following constraint^: 



d^xdy V-G< 5 > (T£ - 2T|) = 0. (3.2) 

To show how this condition actually emerges we work out an example of T M N created 
by a 5D scalar field, noting that the same arguments can be extended for other types of 
matter. We start by making a standard Kaluza-Klein decomposition of our diagonal metric 
( |2.2j ) in the following way 

G MN = ^ ^ g 2 7 ^ > ( 3 ' 3 ) 

where Gffi denotes the usual 4-dimensional metric and e 2 " 1 ^ = b 2 (t). Here, for the sake of 
our analysis, we restore the time-dependence of the scale factor along the extra dimension. 



1 This constraint was also derived from a topological argument in [pT| 
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We also consider the following 5-dimensional field theory of a single scalar field minimally 
coupled to gravity 



S = -J d 4 xdy^f^l^ + ^d M <t>d M ( f> + 



V(cf>) 



(3.4) 



where k 2 = 87r/M|. Under the aforementioned decomposition ( |3.3|) , the above theory can 
be written as 



S 



J d 4 x dy V^GW e< | ^ [r® (t,y) - 2e^ D^e? 



1 



+ ^d M <pd M <P + V(<P)\ (3.5) 



The above form of the theory is not particularly convenient for our analysis: the "kinetic" 
term for the dilaton field 7 is a total derivative, and, thus, drops out from the action and, 
moreover, there is a direct coupling between 7 and R^(t,y). Both of these problems can 
be overcome by making the conformal transformation Gffi = e~ 1 g^ iU . Then, the action 
functional of the theory takes the final form 



S = -J d i xdy^^[R^(t,y) + h fll d^ 



(3.6) 

Our next task is to derive the equation of motion for the dilaton field 7 and find the 
condition that the components of the energy-momentum tensor of matter need to satisfy in 
order to ensure the stability of the extra dimension. The variation of the action ( |3.6|) with 
respect to 7 is given by the expression 



5S 



-37 (^f- e ^)]U 7l 



from which we learn that the stability condition, analogous to (|3.1|) is 



dyV^[-e-^ (d 5 0) 2 + e-n^(0) 



0. 



(3.7) 



(3.f 



This constraint can be expressed in terms of the components of the energy-momentum 
tensor generated by the 5-dimensional scalar field. Before making the conformal transfor- 
mation of the 4-dimensional metric, Tmn had the form 



T) 



MN 



d M 4> d N (b - G 



(5) 
MN 



d P (t)d p (t) + V((t)) 



(3.9) 



Then, we can easily find the following components 

= G (4)^^ = _ 2e - 2 ^ («9 5 0) 2 - W{<f>) 



e- 27 f 55 



\d,<pd^+ l -e~^ (9 5 0) 2 -V(0). 



(3.10) 
(3.11) 



We can now rewrite the constraint equation ( |3.8|) preventing a linear source for the dilaton 
7 in terms of Tmn which up to a multiplicative constant yields 



e 



-7 



dyV^g — (r;-2T 5 5 ) = o. (3.12) 

Note that, although and T| are not conformally invariant, their combination T^ — 2T| is, 
since the difference does not depend explicitly on Gffi. Therefore, the above result remains 
invariant under the conformal transformation that we impose on the 4-dimensional part of 
the metric tensor. 

We can easily demonstrate that the introduction of any type of matter in the context 
of the 5D theory serves as a source term for the 4D gravitons. More specifically, consider 
the linear term 

5G$ N f MN = 5G$ f ^ + die 2 "') f 55 (3.13) 
in any system. Using G$ = e -7 ^, one can rewrite the above expression as 

5G% ] N f MN = e^5g, u f - 5 1 (T£ - 2T 5 5 ) (3.14) 

Thus, we see that T^ v is a source for the graviton while the combination (T^ — 2T|) is the 
source for the dilaton or for any Brans-Dicke scalar field. In the presence of a mechanism 
that stabilizes the extra dimension, the source for the dilaton is identically zero. If we do not 
have a stabilization mechanism but, nevertheless, we look for solutions with a Newtonian 
limit, we still have to require the absence of the above term. On the other hand, the 
presence of a scalar field in the theory will give rise, not to a standard Einstein theory of 
gravity, but to a scalar-tensor theory with all its consequences. 

All of these considerations help to elucidate the required form for T M N in the 5-dimensio- 
nal space which allows us to keep the size of the extra dimension fixed. Eq. (|3.2|) is the 
constraint which in fact must be satisfied. This constraint can, for example, be easily 
specialized to the cases already considered in previous publications. Following |[L0|| , one 
may choose two branes with opposite matter densities and vanishing T|. In this case the 
stability condition ( |3.2| ) is trivially satisfied (since the T^'s of each brane cancel each other). 
As another illuminating example, we consider the thin-wall cosmological solution that was 



derived in Ref. [|19[] and describes the 5-dimensional space-time around an infinitely thin 
brane-universe located at y — 0. The ordinary matter density is confined to the brane and 
no other form of energy - apart from a smooth distribution of T| - exists in the bulk. The 
energy-momentum tensor of the brane-universe, in this case, has the form 

f M N = diag(^M(-p,p,p,p),f 5 5 ) . (3.15) 

The expression for T| in the bulk follows from the assumption that stable extrema of n(y) 
and a(y) exist outside the brane, at y — ±\y min \, and that the Friedmann equation has the 



8 



usual form at the origin. The resulting energy-momentum tensor component is 

f5 _ qg(t) (-p + 3p) , 

' 5 ~ 2n(t,y)a*(t,y) 2b \y mm \ + U[P > " ^ lbJ 

Note that, in our expansion in powers of p, the prefactor (ag/na 3 ) = 1 + 0(p 2 ), so that 
f| = (-p + 3p)/(2b \y mm \) + 0(p 2 ). 

Going back to eq. (|3.15|) , the 4- dimensional trace of the energy-momentum tensor is 
simply given by the expression 

<%) 
bo 

Taking into account the fact that the 4-dimensional part of the energy-momentum tensor is 
localized on the brane while T| is not and performing the integration in ( [3.12 ) with respect 



f; = -^(-p + 3p). (3.17) 



to y G (— |y m j n |, \y m in\)i we can easily demonstrate that the desirable form of T| satisfies 
our general constraint (|3.2j ). 

The general condition on the stability of the radius Q3.2|) represents a highly nontrivial 
fine-tuning problem for the matter fields in 5D space if one insists on a flat potential for 
the dilaton and on the compactification at an arbitrary value of 7. Moreover, the dilaton 
in this scenario will remain massless, which constitutes another serious phenomenological 
problem as noted above. Therefore, we must consider a more realistic situation with an 
additional physical input, namely an explicit mechanism which gives the mass to the dilaton 
and stabilizes it around a specific local minimum. 

It is tantalizing to write this potential in the simplest possible form, 

V{C$) = a(Gg - bl)\ (3.18) 

so that (755 = e 27 = 6q represents an equilibrium point. Since this form explicitly breaks 
general covariance, it must be only effective, realized by some bulk mechanism, unspecified 
at this point. Any uniform deviation of G55 from 6q induces a uniform value for T| in 

(5) 

the bulk, given by the derivative of V with respect to G\§ . If a brane with the non- 

(5) 

vanishing matter energy density is included, the equilibrium position of G55 changes. We 



next show that the value of T|, generated this way, does satisfy eq. (|3.16Q independent of 



the parameters which characterize the stiffness of the dilaton potential. 

To have a workable framework which preserves general covariance we choose an auxiliary 
bulk scalar \ with the following unusual action 

rf 4 x f lym " A JI&V dy (d MX d M X - elf (3.19) 

" I Vmin I 

in the given background configuration x{v) = c 2\y\- This background configuration creates 
the potential of the form (|3.18|) if we identify b = c 2 /ci. The dilaton chooses its preferred 
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point 70 and it is easy to check that G55 = e 2 " 70 = feg = c 2/ c i together with x{y) — c 2\y\ 
satisfy all of the equations of motions. At the same time, at the equilibrium position, all the 
components of T M N which correspond to the action (|3.19|) remain identically zero. We also 



note here that such a background configuration for x is consistent with compactification 



and it can be further checked that with the choice of (|3.19|) , the cusps in x(u) do not require 
the existence of delta-functional sources. 

When a brane with matter content is introduced at y — 0, the equilibrium position for 
7 changes. This change can be read off from the interaction of 7 with the trace of the 
brane energy- momentum tensor, eq. ( |3.12| ), and the effective potential generated by ( |3.19|) . 



Keeping only the term linear in this deviation, 7 — 70, in front of the small perturbation 
created by of the brane, we write down an effective potential for the dilaton of the form 

8 1 

Veffd) = \ymin\ 4(1 - 7o) 2 - (-P + 3 P)(l ~ 7b), (3-20) 

from which we see that the new equilibrium point is given by 

(7 - 70) = -r o 2 [7 I • ( 3 - 21 ) 

c 2 OZj \ymin\ 

This deviation produces a non-vanishing (55)-component of the energy-momentum tensor 
in the bulk (from the variation of |3.19| ), exactly of the form discussed above. Indeed, 



neglecting all terms on the order of p 2 , we have the following value for T| 

f 5 5 = 84e"^(7 - 70) = 7TT^7> (3-22) 

whereas the (00) and (ii) components of T M N in the bulk remain zero to this accuracy. The 
same arguments hold for any mechanism which would render an effective potential in the 
form ( grp . 

Equation ( 3.22| ), is exactly of the same value needed for the restoration of the Friedmann 
equation for the brane scale factor, and allows for an important interpretation of T|. Indeed, 
we saw that this component arises as the backreaction of the dilaton fixing potential on 
the presence of the energy density on the brane. This confirms a hypothesis put forward in 
our previous work [0. At the same time it removes the necessity of the fine tuning (|3.12| ) 



required for the case of the flat potential. It turns out that the only requirements which 
should be imposed on a in ( |3 . 1 8|) or c 2 , in (|3.19|) - (|3.22|) is the condition 7 — 70 70 which 



justifies our linearized approach. We see that the details of the dilaton fixing potential are 
irrelevant to linear order in p and affect only the sub-leading, p 2 -proportional terms. We 
note that mechanisms for generating a dilaton potential from bulk fields was discussed in 
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Recent work on the stabilization of the dilaton in the context of the Randall- Sundrum 
model showed a similar shift in the dilaton expectation value |p3f . Indeed, the shift found 
there and the resulting change in the derivative of the dilaton potential at the shifted point 



acts as a source to G55 giving an identical expression for T 55 described here and in [19 



4 Solutions with cosmological constant and compact- 
ification 

In this section, we adopt the scenario of an infinitely thin brane-universe located at y = 
with the ordinary energy density and pressure being localized on the brane. We further 
assume that the total energy density of the brane, pi, r , is the sum of the ordinary matter 
density, p, and a brane cosmological constant, A^. On the other hand, T| is smoothly 
distributed along the extra dimension while the bulk is dominated by the energy density 
coming from a bulk cosmological constant, A#. Here, we are going to investigate the ex- 
istence of cosmological solutions that admit compactification of the extra dimension while 
leading, once again, to the standard cosmological expansion of the brane scale factor. We 
will consider three different cases corresponding to negative, positive and zero bulk cosmo- 
logical constant. The sign of A{, r will be determined by demanding the cancellation of the 
effective cosmological constant in the Friedmann equation. Related brane solutions with 



both brane and bulk cosmological constant were considered in [JT2|, [24 |. 



4.1 Negative cosmological constant in the bulk 

We start with the case of a negative bulk cosmological constant, A# < 0, and we write the 
energy-momentum tensor in the bulk in the following form 

T£ = (-Pb,Pb,Pb,Pb,Ts) , (4.1) 

where p# = A# < and the other components are, for the moment, arbitrary. For a 
smooth, constant distribution of energy in the bulk, the solution for the scale factor a(t, y), 
outside the brane, will follow from the (OO)-component of Einstein's equations, which takes 
the simple form 

x" - A 2 x - B 2 = , (4.2) 

where x = a 2 and 



The general solution of the above differential equation has the form 

A 2 



a 2 (t,y) = d 1 (t) cosh(AM) +d 2 (t) smh(A\y\) - ^ . (4.4) 
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The boundary conditions on the brane will determine the two unknown functions of time, 
di(t) and d,2{t). Defining a 2 (t,y = 0) = ag(t), we can write 

di(t) = al{t) + ^ (4.5) 
while the jump in the first derivative of a(t, y) across the origin, which now takes the form 

\n>] k 2 

-— (p + Abr), (4.6) 



a b 3 

leads to the following expression for the second function of time 



d 2 (t) = (p + A br ) . (4.7) 



When the general solution ( |4.4| ) is substituted in the (n)-component of Einstein's equa- 
tions, we get the anticipated result p B = |A_b| = — ps while, the same procedure with the 
(55)-component, leads to the constraint 



d_ 

dt 



A 2 , , Q , , b A 



2 



^6Vd(f|-|A B |). (4.8) 



The above constraint will allow us to determine the form of T| once the expression for X(t) 
is fixed by the Friedmann equation. 

Next, we come to the question of the compactification. Even in the absence of a second 
brane, the compactification of the extra dimension could easily take place as long as the 
solution (|4.4j ) has a stable extremum at some point y = \y m in\- Then, by identifying the 
points y = ±\y m i n \, the extra dimension would be compactified with its size being 2b \y m in\- 
From the vanishing of a'(t,y), we find that 

tanh(A \y min \) = — ~7~7~r = c . (4.9) 
ai{t) 

The above ratio, as indicated, must be constant in order for the extremum to be stable in 
time. Then, by using the above relation and eqs. (|4.5|) and (|4.7|) , the Friedmann equation 
immediately follows 

1 = ^|A B |{_ 1 + ^i^±M}.^. (4.10) 
oo/ A 2 ag 3 I c JQk 2 \A b \ ' 3 



The last equality in the above equation defines the 4D Planck constant and holds only 
under the constraint 

_ i + / L A^ = Q (4n) 



c o J6k 2 \A B \ 
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and the redefinition 



2 k\lki\A B \ 



c V 6 



These expressions can also be used to define the 4D cosmological constant on the brane. 
From eq. (|4.10 ), we can write 



co J6k 2 \A b \ 



-V -'W^r ( 4 - 13 ) 



In the limit that -A|y m i n | is large, Co — 1 and A e // = A& r — y 6 ^ 2 a L In the other limiting 
case, when A|y m j n | is small, we have A e jj = A fer — (26 |y m i„,|)A B . 

From its definition ( |4.9| ), since A|y m i n | is positive, it is easy to see that cq is always 
positive independently of the exact values of |A_b| and \y m in\ and, thus, the Friedmann 
equation will always have the correct sign. Then, in order for the constraint (|4.11|) to be 
satisfied, we must necessarily have A^ > 0. The above two equations can be rewritten in 
a simple form as 

|ABl = 6^ A ' r ' K2 = i? Abr - (414) 

Finally, we need to determine the form of T| in the bulk. By using the constraint 
and the relations ([4.14 ), we find that 



A number of comments should be made at this point. In the limit \y m i n \ — > 00, the extra 
dimension is non-compact and Cq — 1. In that case, the last term in the above expression 
becomes zero and, modulo a small correction of 0(p 2 ), T| assumes the form of another 
pressure-like component of the energy-momentum tensor, equal to \Ag\. This is the case 
considered in Ref. ||24|| . In order to compactify, their solution requires the introduction of 



a second brane along with an unavoidable correlation between the energy densities on the 
two branes. In particular, the absence of T|, again implies the presence of a negative energy 
brane. For every other value of cq, our results describe 5-dimensional cosmological solutions 
that respect the 4-dimensional Friedmann equation and allow for the compactification of the 
fifth dimension via the existence of an extremum of a(t,y) at some finite point y = \y m %n\- 
As we have seen in the previous section, the deviation of the fifth component of the energy- 
momentum tensor from the strict limits of a pressure-like component, is crucial for the 
stabilization of the extra dimension. 
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Let us, finally, note that the compactification of the extra dimension via the extremum 
of the scale factor not only eliminates the need for the introduction of a second brane but 
also resolves a major problem that arises in the framework of the two-brane models. For 
Ab < 0, it has been shown |13|, 14] that, in order to reproduce the correct sign in the 
Friedmann equation, one of the two branes must have an energy density which is negative, 
p < 0. The need for this unnatural assumption ceases to exist in the context of our analysis. 
Despite the fact that we start with the same assumption of a negative cosmological constant 
in the bulk, in our case, the correct sign in the Friedmann equation is always maintained. 
Two- and multi-brane solutions in this context will be presented elsewhere. 



4.2 Positive cosmological constant in the bulk 

In this case, we assume that p# = A# > 0. The analysis closely follows the one of the 
previous subsection, however, the opposite sign of A B modifies some of our results and, 
consequently, our final conclusions. We start again with the general solution for the scale 
factor a(t, y) outside the brane, which, now, has the form 

B 2 

a 2 {t, y) = cos{A\y\) + d 2 (t) sm{A\y\) + . (4.16) 

The first coefficient d\{t) is given by eq. ([4.5|) with A 2 being replaced by —A 2 while the 
second one, <i 2 (t), is still given by eq. ( |4.7| ). Once again, it is easy to see that the above 
solution satisfies the (n)-component of Einstein's equations provided that ps = —Ab = ~Pb 
while, the (55)-component, leads to the constraint 



d 
dl 



A 



2 



(4 + d( 



4 v ~ z ' ~ LJ A 2 \ A 



2k 



2 



6Va(T 5 5 + A B ). (4.17) 



We still need a stable extremum for the compactification of the extra dimension and 
this extremum, now, takes place at 

ttm(A \y min \) = = c . (4.18) 

By using the above relation and the expressions for di(t), the Friedmann equation can, now, 
be written as 

'd \ 2 k 2 f k 2 (p + Abr) 



a j 3 I c V6rA B J 
The standard form of the above equation is immediately restored under the constraint 

c v6ft A B 
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while the redefinition of the 4-dimensional Newton's constant in terms of the 5-dimensional 
one is still given by eq. ( 4.12|) . From its definition, eq. (|4.18 ), we conclude that Co is not 



always positive. Thus, in order to have the correct sign in the Friedmann equation, we have 
to impose the following constraint on the possible values of \y m in\ 

U7T < A\y min \ < (n + -)7T, neZ. (4.21) 



Then, the constraint ( |4.20|) can be satisfied only in the case where Ab r < 0. Note that the 



cosmological constants in the bulk and on the brane must always have opposite signs in 
order to ensure the cancellation of their contributions to the energy density of the universe. 
However, their sign does not enter the Friedmann equation neither imposes any unnatural 
constraints on the energy density of our brane-uni verse. 

Let us note that, contrary to the case where A# < 0, here, there are no solutions 
that would correspond to a non-compact extra dimension. Every cosmological solution, 
described by eqs. ( [4.1 6| ) and ( |4.18| ), that respects the Friedmann equation at the origin 



has a compact extra dimension with finite size 2b\y m i n \. This is also obvious from the 
expression for the fifth component of the energy-momentum tensor in the bulk, which from 
the constraint (14.171 ) is found to be 



= " A * + 12n(S (< ,y) { " ' ( " + 3P) + + 3) A " <2Air + ' ~ '"'I ' (4 ' 22) 

Note that the form of T| is always different from the one of a pressure-like component and 
there is no real value of Co that would eliminate the last term. Thus, in order to stabilize 
the extra, compact dimension, a non-trivial T| must be generated. 



4.3 Zero cosmological constant in the bulk 

Finally, we consider the case where A# = while keeping a non-zero cosmological constant 
on the brane. In the bulk, the solution is exactly the same as in the case studied in the 
thin-wall approximation of our previous work |]TP[. The general solution for a(t,y) can be 
written as 

a 2 (t,y) = a 2 (t) + c(t)\y\ + ^y\ (4.23) 

with 

k 2 

c(t) = -—a 2 b(p + A br ). (4.24) 

For p + Ab r > 0, the above solution describes a spatial scale factor that decreases as we 
move away from the brane. However, the rate of decrease is much smaller than the one 
given by the warp factor in the Randall- Sundrum model 0. 
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The compactification of the extra dimension is made by identifying the two extrema 
which, now, take place at the value 



_ dtm) 

\ymin\ — ' V ± ' Z " J J 

and, by using the above result for c(t), the Friedmann equation may be written as 

+ ^(p + A 6r ). (4.26) 



a / 3 2b\y min \ 3 

Due to the absence of a bulk cosmological constant that would cancel the contribution of 
A&r to the total energy of the universe, the term linear in p will become dominant only in 
the case Ab r <C p. In the opposite case, as anticipated, the presence of Ab r will significantly 
modify the expansion rate of the scale factor at the origin. 



5 Thick-wall solution with compactification 

To complete the class of solutions we have been discussing, we generalize our thin-wall 



solution found in [19j to one of finite thickness. We assume that our brane-universe has a 
non-zero thickness 2A6 with the normal matter content being distributed throughout the 
brane. Outside the brane, both the energy density and pressure are zero, which leads, after 
integrating the (OO)-component of Einstein's equations, to the following general solution for 
the spatial scale factor 

a 2 (t, y) = al(t) + C(t) (\y\ - A) + ^ (y 2 - A 2 ) , (5.1) 

where a&(t) is the scale factor evaluated at y = A and C(t) an unknown function of time 
which needs to be determined. For the compactification of the extra dimension, we require 
an extremum in a(t,y) outside the brane, so, we impose the condition a'(t,y min ) = on 
the general solution (|5.1|) . That leads to the result 

(52) 

\ymin\ ' V u "^7 



In Ref. |[I9|| , we showed that, in the thin- wall approximation, the existence of a stable 
minimum for a(t, y) outside the brane and, thus, the compactification of the extra dimension 
is equivalent to the restoration of the Friedmann equation at the origin. Here, we need to 
impose that this equation takes the usual form, i.e. 

w = i*m (5 . 3) 

a J \ 2 a,Q 3 
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where p(t, 0) is the usual 4- dimensional energy density defined at the origin, and 
Then, by combining eqs. (|5.2|) and (|5.3|), we obtain 



C(t) = -^ (t '° )fl ° & . (5.5) 
3 

It is easy to see that the solution ( |5.1| ) trivially satisfies the (ii)-component of Einstein's 
equations while, when substituted in the (55)-component, it leads to the constraint 

d T b 2 , A , ls o cia 1 2/t 2 o. , . 

= a 6 aTg . (5.6) 



dt 



^4 IZ/minl) ^ 2 



The above constraint, with the expression for A given by eq. ( |5.3|) , is satisfied only if the 
fifth component of the energy-momentum tensor, in the bulk, has the form 

"It (B) (* , y) = A { 7^7 ["A (p+3p) - 2 P a n A ) - ^ (A ~ ^ ? P (P+3p) } , (5-7) 



where again p and p are defined at the origin. In the limit A — > 0, the above result reduces 
to what we had for the thin- wall solution [[19]. We emphasize once more, that despite the 
peculiar nature of the expression in (|5.7|) , so long as the dilaton picks up a stable expectation 
value, this is the form that T 55 will assume. 

To complete the thick wall solution, we must derive the expression for T 5 5 (t,y) inside 
the brane. This derivation is given in the Appendix, the result is 



Ti {br) (t,y)=fi {B) (t,y)+p 



n(t, A) 
n(t,y) 



(5.8) 



where T^r B Jt,y) is given by eq. (|5.7|). We can easily show that the above expression 
satisfies the fifth component of the equation for the conservation of energy ( |2.9p and that, 
in the limit y — > A, the second term at the r.h.s. of eq. (|5.8|) vanishes, thus, ensuring the 
continuity of T| across the boundary. 

There is one remaining function we must determine, and that is the scale factor on the 
boundary of the brane, a A (t). Again, we leave the derivation to Appendix, and simply 
display the result here 

a A (t) = ao(t) [l - — p(t,0)Ab\ . (5.9) 

The above result reveals the fact that the scale factor a(t, y) decreases as we move away 
from the origin. Nevertheless, the rate of decrease being simply proportional to the energy 
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density on the brane is much smaller than in the case of the Randall- Sundrum thin- wall 
solution || where the rate of decrease depends exponentially on the bulk cosmological 
constant. Due to the continuity condition on a'(t,y) at the boundary, the scale factor will 
also decrease mildly outside the brane until it reaches its unique minimum at y = ±\y min \. 
The mild decrease of the scale factor makes necessary the compactification of the extra 
dimension which is realized through the identification of the two stable extrema ±|y m j„|. 

In this section, we have demonstrated that the positive features exhibited by the thin- 
wall solution found in |19| also characterize its thick-wall analog. The cosmological solutions 
that emerge respect the 4-dimensional Friedmann equation on the brane and the compact- 
ification of the extra dimension takes place without the need for a second brane. In both 
instances, the existence of a non-vanishing Tj plays an important role as it results from the 
stabilization the extra dimension and helps restore standard cosmology on the brane. 



6 Conclusions 

While it may seem that standard FRW cosmology should arise trivially on a homogeneous 
3-brane embedded in a higher dimensional space time, it has become clear that this result is 
highly dependent on the generation of a dilaton potential which fixes the dilaton expectation 
value and hence the size of the extra dimension. In the absence of a dilaton potential, there 
are two serious problems which must be overcome. Aside from the non-compact solution 
of Ref. ||, compactification of the extra-dimension in solutions for 3-branes with matter 
require the presence of a second brane on which the matter is constrained by Einstein's 
equations to take the form p 2 = — Pi + 0(p 2 ). That is, we are forced to tolerate negative 
energy densities on the second brane. This is most assuredly unphysical. However, even 
if we accept this solution which includes the negative-brane, the solutions to Einstein's 
equations yield a cosmological expansion for which the Hubble parameter is proportional 
to the energy density. BBN and post BBN-cosmology can not be reconciled with this 
expansion law. We stress that we distinguish here the difference between a brane with 
negative tension and one with negative matter density. While a negative tension brane 
may naturally appear in string theory, we are not aware of any formulation which allows 
for a negative matter density —p. 

In our previous work, we proposed that both of the aforementioned problems were tied 
to explicit assumptions made on the form of the energy-momentum tensor. Namely, that 

f M N = 6 -P-dmg(-p, P ,p,p,0) . (6.10) 

i.e., that the (55)-component of the energy- momentum tensor vanished. In |JL9[| , we showed 
that by allowing T55 to differ from zero, we could find a solution for the 3-space scale factor 
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a(t,y), in the bulk, which has a local minimum at some co-ordinate value y m in, at which 
point we can compactify the extra dimension without the introduction of a second brane 
and in particular without a negative energy brane of any kind. Secondly, we derived the 
form of T|, in order to recover a normal Friedmann expansion on the 3-brane. 

Here we have justified the choice for T|, and showed that it is the natural result of 
the back-reaction of a stabilized dilaton to the 3-brane with matter. In particular we 
have showed that the dilaton couples to the combination — 2T|, and any stable dilaton 
configuration requires this combination of energy-momentum components to vanish when 
averaged over the extra dimension. Furthermore, when a dilaton potential, which fixes the 
dilaton vev is included in the derivation of the solution, matter on the 3-brane will induce a 
shift in the dilaton vev in such a way so as to exactly guarantee the vanishing of — 2T|, 
and produce the needed (55)-component derived in [[0|. 

We have, in addition, supplied generalization of our 3-brane solutions to include a cos- 
mological constant both in the bulk and on the brane. We have considered the cases of 
negative, positive and zero bulk cosmological constant and demonstrated that a solution 
with a stable, compact extra dimension, which respects the standard 4D Friedmann equa- 
tion, arises in each case. Moreover, the correct sign in the Friedmann equation is always 
maintained for every sign of as long as the sign of A& r is exactly opposite. In this case, 
the contributions of the two cosmological constants to the energy density of the universe 
can be cancelled leaving a term linear in p to govern the expansion of the brane scale fac- 
tor. Finally, our single-brane thin-wall solution derived in [HJ was extended in the case 
of a brane-universe with a finite thickness. The resulting cosmological solution allows for 
the compactification of the extra dimension through the identification of the two stable 
minima of the scale factor outside the brane and the 4D Friedmann equation is once again 
recovered. In all of the above cases, the existence of a non-trivial fifth component of the 
energy-momentum tensor is necessary for the stabilization of the size of the extra dimension. 
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Appendix 

In this Appendix, we show the details of the derivation of the energy-momentum tensor 
inside the brane and subsequently generalize this solution for an arbitrary equation of state. 
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For reasons, that will become obvious, we may rewrite the expression for the scale factor 
outside the brane (15. If) as follows 



a\t,y) 



A 2 (t) 



where 



C out (t) = —C(t) = -\y r: 



\y\ + C out {t)} 



E out (t) = 2b 2 



2 X 2 (t)E out (t) 



2b 2 



y ( A - \Vr, 



\2 a A 



,2 
A 

A 2 



(A.1) 



(A.2) 



To obtain the inside solution for the scale factor a(t, y), we will assume, at first, that the 
pressure on the brane is zero, so our results would apply for a matter-dominated universe. 
Then, the zeroth-component of the equation for the conservation of energy fl2.8|) gives 
P = Po/ fl3 - Integrating the (OO)-component of Einstein's equations, the general solution for 
the scale factor takes the implicit form 



2A 2 



lo 



B 3 (t) b \B(t) 



B 2 (t)a(t, y) - A 2 \ + 2^ E m (t) + B 2 (t)a 2 (t, y) - 2A 2 a(t, y) 
2 



+ 



B\t) 



E m (t) + B 2 (t)a 2 (t, y) - 2A 2 a(t, y) = ±V2 [\y\ + C m {t)} , (A.3) 



where 



A 



2b 2 k z p 



B 2 (t) 



2b 2 

W)' 



(A.4) 



This is the same expression as the one that we derived in the thick- wall approximation [19 



for the scale factor on the brane. The function Ci n (t) can be determined by evaluating the 
solution at y = 0, and has the form 



Ci n (t) 



V2A 2 



log 



B*(t) °\B(t) 



B 2 (t)a (t) - A 2 



(A.5) 



The remaining unknown function Ei n (t) will be determined from the continuity of the first 
derivative of a(t, y) with respect to y at the boundary y = A between the brane and the 
bulk universe. If we use the form (|A. 1|) for the solution outside the brane, this condition 
can be written as 



± V2a A a' A = yJE^t) + B 2 {t)a\ - 2A 2 a A = y/ ' E out {t) + B 2 {t)a\ (A.6) 
leading to the result 



E in (t) = 2b 2 



(A - \y min \f 



A 5 



+ 



2k 2 



Po a a 



(A.7) 
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Moreover, the function E{ n (t) is constrained by the (55)-component of Einstein's equations. 
When the solution (|A.3|) is substituted in eq. (|2.6| ), we obtain the following constraint 



dE in (t) Ab 2 



k 2 a 6 a (T 5 5 + p) 



(A. 



dt 3 

This is satisfied if, and only if, the fifth component of the energy- momentum tensor, T|, on 
the brane is given by the expression in ( |5.8| ). 

We have one more condition that our solution on the brane needs to satisfy, namely, 
the vanishing of a'(t,y) at the origin. This leads to the result 



E in {t) = 2A 2 a (t)-B 2 (t)a 2 l (t) 



(A.9) 



2a \y„ 



A (a + a A ) 



(A.10) 



If we substitute the expression for E in {t) ( |A.7| ) in the above equation, we find the following 
constraint 

— (A — \ 1 , h 2 _ (ao-flA) 

^2 V \ymin\) ^ 

where we have set po = Po/2A5. In the limit A — > 0, the above constraint should reduce 
to the expression (|5.2|) for \y m i n \- This condition allows us to determine the function a&(t) 
given in ( {5.9|) . 

Finally, we extend the above analysis for the determination of the thick-wall solution 
inside the brane, for a matter-dominated universe, to the case of a brane-universe with a 
general equation of state p = w p. In this case, the equation for the conservation of energy 
gives the result p = po/a 3 ^ 1+w ' . The first derivative of the scale factor on the brane, now, 
takes the form 



a\t,y)=± 



lE in {t) + B 2 (t)a 2 (t,y) 



A 2 a(t,2/)( 1 - 3 "0 



V2a V""" ' " y '~ v " a/ (l-3w) 
and its continuity across the boundary y = A leads to the following result 



E m {t) = 2b 2 



r h2 



A 4 



(A -| Vmin 



2k 2 



X 2 



3 (1 - 3w) 



(l-3w) 

p a\ 



(A.11) 



(A.12) 



The constraint ( |A.8| ) that follows from the (55)-component of Einstein's equations still 
holds and leads to the following expression for the value of T| on the brane 



fl {br) (t,y)=fi {B) (t,y)+p 



n(t,A) ( a(t, y y 
n{t,y) Va(^ A ), 



3w 



(A.13) 



As in the previous case, we still have to fulfill the vanishing condition of a'(t,y) at the 
origin. Then, we find the result 



b , » i >o 1 J 2 1 y m in | G^o 



(A -I ymin | / 



A 



3w 



1 _ ( a ° \ 3w ~ 1 



- A (oq - a^) 



(A.14) 
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which, in the limit A — > 0, reduces to 



v a ' 



12 



(A.15) 



Note that for every value of w, larger or smaller that 1/3, the scale factor decreases as we 
move away from the origin, thus, ensuring the monotonic behavior from the point y = to 
its minimum at y — \y m in\- 

However, the above expression gives a trivial result for the choice w = 1/3 which 
corresponds to a radiation-dominated universe. Moreover, the expression ( A.ll|) for the 
first derivative of a(t,y) diverges for the same choice. This means that we have to study 
this case separately. Starting all over again and setting p = p /a 4 , the first derivative, 
a'(t,y), takes the form 



a'(t, y) = ± — JE in (t) + B 2 (t)a 2 {t, y) - 2A* lna(t, y) 
The result for E in (t) is accordingly modified and is given by 



E in (t) = 2b 2 



n 2 ?k 2 

w / A I l\2 A a 1 

—r (A - \y min \) - — + — p lna 



A 



A 2 



(A.16) 



(A.17) 



leading to the following form of T| on the brane 



rp5 

1 5 (br) 



T 



\b)+P 



n(t,A) a(t,y) 



n(t,y) a(t,A) 

Finally, the vanishing condition on a'(t,y) at the origin leads to the constraint 



In the limit A — > 0, we obtain the result 



(A.18) 



^ (A - \y mm \) 2 = i ^2\y min \a 2 In (^-) - A (a 2 - J . (A.19) 



(t) = a (t)exp(-^p(t,0) A6). 



(A.20) 



Once again, the scale factor on the brane is decreasing, as y increases, resulting in a similar 
behavior outside the brane until its minimum. 
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